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Введение. Рассмотрен перспективный способ обеззаражива- 
ния стоков животноводства. Суть способа заключается в ком- 
плексном физико-химическом воздействии переменного маг- 
нитного поля в сочетании с химическим реагентом. Устрой- 
ством реализации предложенного способа является индуктор. 
В качестве параметра, характеризующего энергетическую 
эффективность технологического процесса обеззараживания 
стоков животноводства, выбрана удельная электрическая 
энергоёмкость. Для повышения эффективности исследования 
технологического процесса обеззараживания стоков животно- 
водства предлагается использовать математическую теорию 
планирования эксперимента. 

Материалы и методы. В ходе проведения данного исследо- 
вания использовались методы: инструментального, прямого и 
косвенного измерения, статистической обработки данных, 
математическая теория планирования эксперимента. 
Результаты исследования. Экспериментально исследовано 
влияние технологических факторов на энергетическую эф- 
фективность технологического процесса обеззараживания 
стоков животноводства. На основании полученных данных 
выполнено ранжирование по силе влияния технологических 
факторов на удельную электрическую энергоёмкость процес- 
са. Установлена несущественная степень влияния оценок ко- 
эффициентов парного взаимодействия, на что указывает их 
отсутствие в математической модели. 

Обсуждение и заключения. Создана математическая модель 
технологического процесса обеззараживания стоков животно- 
водства, позволяющая оценить влияния технологических фак- 
торов на удельную электрическую энергоёмкость процесса 
обеззараживания. 


Ключевые слова: обеззараживание, стоки животноводства, 
индуктор, математическая теория планирования эксперимен- 
та, математическая модель. 
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Введение. Одним ИЗ ОСНОВНЫХ факторов загрязнения окружающей среды являются промышленные стоки 
сельского хозяйства. Наибольшую опасность представляют сточные воды животноводства (далее стоки), являющиеся 
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источником большого количества заболеваний, утилизация которых невозможна без их обеззараживания. Наиболее 
перспективным способом обеззараживания стоков является комплексное физико-химическое воздействие вращающе- 
гося переменного магнитного поля и активного хлора [1]. Энергетическим параметром качества процесса обеззаражи- 
вания является удельная энергоёмкость процесса. Соответственно, создание математической модели, позволяющей 
оценить влияние технологических факторов на энергоёмкость процесса, является актуальной задачей. 

Решить эту проблему можно используя математическую теорию планирования эксперимента. 

Целью данного исследования являлось создание математической модели технологического процесса обеззаражи- 
вания, оценивающей его энергетическую эффективность. 

Экспериментальное исследование. Для оценки энергетической эффективности процесса необходимо рас- 
сматривать затраты электрической энергии на получение единицы готовой продукции, то есть удельной электриче- 
ской энергии. В качестве единицы готовой продукции лучше всего рассматривать её объём. 

На основании вышеизложенного, при оценке энергии, затраченной на технологический процесс, необходимо 
рассматривать активную составляющую удельной электрической энергии, расходуемой на получение единицы гото- 
вой продукции. 

Представленные в табл. 1 факторы были приняты в качестве способных оказывать влияние на технологический 
процесс обеззараживания стоков комплексным физико-химическим воздействием в индукторе. 























Таблица 1 
ТаЫе 1 
Факторы и уровни их варьирования в натуральных переменных 
Каслтогз апа [еуе5 ор йе" уатаБИйу т рйуяса! уапа е5 
Кодовое Уровни факторов 
Фактор 
обозначение Х=- Х=О | ХЕ 
Заполненность рабочими телами рабочей зоны индуктора р в % Х! 0,74 2,96 5,18 
Отношение длины рабочих тел к их диаметру /4 А? 2) 15 25 
Магнитная индукция В в мТл ХЗ 40 60 80 
Концентрация активного хлора в в мг/л А 6 12 18 
Продолжительность воздействия {вс Хх 2 4 6 




















В качестве условий исследования рассматривались [2—5]: 

— частота колебаний электрического тока в сети /, Гц; 

— материал рабочих тел. 

Кибернетическая модель системы (факторы-условия-объект исследования-параметры показатели качества), ил- 


люстрирующая их взаимосвязь, представлена на рис. 1. 
Факторы Условия Показатели качества 


Й Материал 
ь рабочих тел 
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4 Технологический процесс 
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2 Рис. 1. Схема кибернетической модели системы 

В при исследовании влияния параметров технологического процесса 

Е обеззараживания стоков на его качество: 
ее р — заполненность рабочими телами рабочей зоны индуктора; //4 — отношение длины рабочих тел к их диаметру; В — магнитная 
2 индукция; © — концентрация активного хлора; #— продолжительность совместного воздействия магнитного поля, активного хлора 
8 и рабочих тел; Г — частота колебаний электрического тока в сети; материал рабочих тел; 

— №>› — удельная электрическая энергобмкость процесса 
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Экспериментальные исследования проводились на индукторе, устройство и принцип действия которого описаны 
в [2]. 

Анализ ИСТОЧНИКОВ [6, 7] и практика решения задач с помощью методов планирования эксперимента показыва- 
ЮТ, ЧТО В большинстве случаев процесс обеззараживания невозможно описать математической моделью первого по- 
рядка ввиду ее неадекватности. Значительная кривизна поверхности отклика, как правило, не позволяет описывать её 
уравнением первого порядка даже в весьма узкой области факторного пространства. В таких случаях необходимо 


описывать исследуемый процесс с помощью полиномов более высоких порядков, например, второго, имеющего вид: 


[ [ к 
| у 2 . 

У=Ь-+ хр + Вим + Рух 
1= 1= = 





где Бо, Бь Би, Б; — коэффициенты, характеризующие силу влияния свободных, 
линейных, квадратичных и эффектов взаимодействия членов уравнения. 

Существует большое количество планов эксперимента, решающих определённые задачи. При выборе вида плана 
для конкретного эксперимента необходимо использовать характерный показатель решения задачи, по значению кото- 
рого оценивается оптимальность найденного решения, то есть, максимальное удовлетворение поставленных требова- 
ний. Таким показателем является критерий оптимальности плана — мера эффективности плана. 

В настоящее время используется свыше 20 различных критериев оптимальности планов [8], которые подразде- 
ляются на две основные группы. К первой группе относят критерии, связанные с ошибками оценок коэффициентов, ко 
второй относятся критерии, связанные с ошибкой оценки поверхности отклика. 

Ввиду того, что целью следующего этапа исследования является оптимизация параметров технологического про- 
цесса обеззараживания, что требует наиболее точное описания поверхности отклика, предпочтительным являются 
критерий С-оптимальности. 

Из всего множества С-оптимальных планов или близких к ним по свойствам, с учётом количества рассматривае- 
мых факторов, наиболее предпочтительным является план Хартли Наз. Этот план имеет минимальное количество 
опытов и предпочтительные статистические характеристики, по сравнению с аналогичными планами, поэтому реко- 
мендуется для создания моделей в виде полинома второго порядка [9, 10]. 

Матрица плана Наз в кодированных переменных представлена в табл. 2. 

Таблица 2 
ТаЫе 2 
Матрица плана и результаты эксперимента 
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Удельная электрическая 
Номер УРОВНЯ ФАЕТОров энергоёмкость 
опыта = ы 

м | № | № | МХ | № Ууээ Ууэ 

1 -1 -1 -1 -1 +1 7,00 7,02 
2 +1 -1 -1 -1 -1 0,20 0,19 
3 -1 +1 -1 -1 -1 3,60 3,59 
4 +1 +1 -1 -1 +1 1,60 1,59 
5 -1 -1 +1 -1 -1 9,60 9,62 
6 +1 -1 +1 -1 +1 7,60 7,62 
7 -1 +1 +1 -1 +1 11,00 11,02 
8 О О ет О ОЙ БС 4.20 4.19 
9 -1 -1 -1 +1 -1 4,60 4,59 
10 +1 -1 -1 1 1 2,60 2,59 
И -1 +1 -1 1 1 6,00 5,99 
12 +1 +1 -1 +1 -1 0,20 0,17 
13 о ие ОЙ Ее пе 12.00 12.02 
14 +1 -1 +1 +1 -1 5,20 5,19 
15 -1 +1 +1 +1 -1 8,60 8,60 
16 1 1 1 1 | 6,60 6,59 
17 -1 0 0 0 0 9,30 9,25 
18 +1 0 0 0 0 4,90 4,96 
19 0 -1 0 0 0 4,60 4,55 
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20 0 +1 0 0 0 3,60 3,66 
21 0 0 -1 0 0 3,10 3,16 
22 0 0 +1 0 0 8,10 8,05 
23 0 0 0 -1 0 5,60 5,55 
24 0 0 0 +1 0 5,60 5,66 
25 0 0 0 0 -1 4,40 4,46 
26 0 0 0 0 +1 6,80 6,75 
27 0 0 0 0 0 5,60 5,59 











Согласно принятому плану было реализовано 27 опытов. Для уменьшения систематических ошибок опыты про- 
водились в случайном порядке. Каждый опыт дублировался три раза. 


Полученные в результате опытов средние значения контролируемых параметров представлены в табл. 2. 

Проверка полученных значений оценок дисперсий параллельных опытов на однородность по критерию Кохрена 
показала, что при выбранном уровне значимости © = 0,05 и степенях свободы числителя {=2 и знаменателя 2 
=27 расчетные значения С„ меньше критического значения С, = 


= 0,21 [11]. Следовательно, гипотеза об однородности 
дисперсий параллельных опытов справедлива. 


После определения оценок коэффициентов, их характеристик и исключения функций из уравнения при статисти- 
чески незначимых коэффициентах, уравнение регрессии для удельной электрической энергобмкости принимает вид: 


У,» =5,59—2,14х, —0,44х› +2,24х; +1145 +1,52%2 —1,48х2 (1) 


По полученным уравнениям были рассчитаны значения удельной электрической энергоёмкости } у»», представ- 
ленные в табл. 2. 


Проверка полученных уравнений на адекватность по Е-критерию показала, что при выбранном уровне значимо- 


сти © = 0,05 и соответствующих степенях свободы / и/е расчетные значения А», меньше критического значения А». 
Следовательно, гипотеза об адекватности уравнений справедлива. 


После подстановки переходных коэффициентов и преобразования уравнение регрессии удельной электрической 
энергоёмкости приняло вид: 


^ 


№,» =4,815—2,789р+ 0,400 1/4 +0,1128+0,570:+0,308р? —0,0141/ 4? (2) 


Уравнение (2) даёт возможность рассчитать значения удельной электрической энергоёмкости №, для любых то- 
чек, находящихся внутри области изучаемого факторного пространства. 


На основании экспериментальных данных (табл. 2) были построены графические зависимости, представленные 
на рис. 2и3. 
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Рис. 2. Зависимость удельной электрической энергоёмкости процесса 
обеззараживания от факторов, связанных с рабочими телами 
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Рис. 3. Зависимость удельной электрической энергоёмкости процесса 
обеззараживания от магнитной индукции и продолжительности совместного воздействия магнитного поля и рабочих тел 


Ер. 3. Оерепаепсе о} зресйс @есвтса! епегэу 1еп5ийу ор 
@зтуеспоп ргосез5 оп тагпейс тдаиспНоп апа диганоп офсотЬте4 ехрозиге 
ортагпейс Пе ап4 иотте те@ит 
Анализ результатов. Анализ результатов опытов и расчётов (табл. 2), уравнений (1), (2) и графических зави- 
симостей (рис. 2, 3) показал: 

— уравнение в кодированных (1) и натуральных (2) переменных в виде полинома второй степени адекватно по 
критерию Фишера характеризуют зависимость удельной электрической энергобёмкости от рассматриваемых факторов 
при уровне значимости © = 0,05. Поэтому оно может служить в качестве математической модели процесса обеззара- 
живания стоков; 


При исследовании влияния рассматриваемых факторов в области интереса на удельную электрическую энергоёмкость 
установлено: 
— Статистически не значимыми, в рассматриваемом диапазоне факторов, оказались оценки коэффициентов 


концентрации активного хлора в пробе, на что указывает отсутствие членов уравнения ха (1) или в натуральных пере- 
менных 4 (2); 


— заполненность рабочими телами рабочей зоны индуктора р существенно влияет на удельную электрическую 
энергоёмкость, на что указывают оценки коэффициентов членов уравнения х/ (1) или в натуральных переменных р (2). 
При этом, с возрастанием значений заполненности рабочими телами рабочей зоны, происходит убывание значений 


удельной электрической энергобмкости по квадратичной зависимости, на что указывает существенная оценка коэф- 
г 20) 
фициента х/ (1) или в натуральных переменных р^ (2); 


— отношение длины рабочих тел к их диаметру /[/4 существенно влияет на удельную электрическую энергоём- 
кость, на что указывают оценки коэффициентов членов уравнения х2 (1) или в натуральных переменных //4 (2). При 
этом, с возрастанием значений отношения длины рабочих тел к их диаметру 4, происходит убывание значений 
удельной электрической энергобмкости по квадратичной зависимости, на что указывает существенная оценка коэф- 
фициента а. (1) или в натуральных переменных р? (2); 


— магнитная индукция В существенно влияет на удельную электрическую энергоёмкость, на что указывают 
оценки коэффициентов членов уравнения хз (1) или в натуральных переменных В (2). При этом, с возрастанием значе- 
ний магнитной индукции В, происходит возрастание значений удельной электрической энергоёмкости по линейной 


зависимости, на что указывает несущественная оценка коэффициента ь (1) или в натуральных переменных В? (2); 


— концентрация активного хлора в пробе «не существенно влияет на удельную электрическую энергоёмкость, 
на что указывает отсутствие членов уравнения х4 (1) или в натуральных переменных 0 (2); 

— продолжительность совместного воздействия магнитного поля, активного хлора и рабочих тел на пробу { 
существенно влияет на удельную электрическую энергоёмкость, на что указывают оценки коэффициентов членов 
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уравнения х; (1) или в натуральных переменных ! (2). При этом, с возрастанием значений продолжительности сов- 
местного воздействия магнитного поля, активного хлора и рабочих тел на пробу #, происходит возрастание значений 
удельной электрической энергоёмкости по линейной зависимости, на что указывает несущественная оценка коэффи- 


циента ел (1) или в натуральных переменных # (2); 


— несущественными оказались все оценки коэффициентов парного взаимодействия, что указывает на отсут- 
ствие эффекта парного взаимодействия рассмотренных факторов; 

— сила влияния факторов распределилась следующим образом: магнитная индукция; заполненность рабочими 
телами рабочей зоны индуктора; продолжительность совместного воздействия магнитного поля, активного хлора и 
рабочих тел на пробу; отношение длины рабочих тел к их диаметру. 

Заключение. Анализ результатов опытов (табл. 2), расчётов (табл. 2), уравнений (1), (2) и графических зави- 
симостей (рис. 2, 3) позволили оценить влияние ряда технологических факторов на удельную электрическую энерго- 
ёмкость процесса, а также наметить гипотезы, объясняющие физическую сущность процессов, происходящих при 
обеззараживании стоков животноводства путём комплексного физико-химического воздействия магнитного поля и 
активного хлора. 
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